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Abstract — A  new  receiver  design  using  iterative  block  decision- 
feedback  equalizer  (BDFE)  is  proposed  for  high  data  rate 
single-carrier  multiple-input,  multiple-output  (MIMO)  underwa¬ 
ter  acoustic  (UWA)  communications.  The  adoption  of  BDFE 
enables  a  sequence-based  log-likelihood  ratio  (LLR)  calculation 
during  iterative  equalization,  thus  leading  to  better  performance 
and  faster  convergence  than  existing  low-complexity  iterative 
equalization  methods  using  symbol-based  LLR  evaluation.  The 
proposed  BDFE  method  is  applied  to  overlapped  blocks  to  reduce 
performance  degradation  at  the  tail  of  each  block.  The  block  size 
is  flexibly  selected  depending  on  the  practical  channel  condition. 
The  proposed  receiver  scheme  has  been  tested  by  extensive  exper¬ 
imental  data  and  proved  to  be  robust  to  different  transmission 
environments  with  consistently  good  detection  performance.  Data 
from  collected  during  the  SPACE08  experiment  near  Martha’s 
Vineyard  and  the  GOMEX08  experiment  at  the  Gulf  of  Mexico, 
are  both  presented  here. 

I.  Introduction 

Underwater  acoustic  (UWA)  communications  are  challeng¬ 
ing  for  three  main  reasons:  first,  the  available  channel  band¬ 
width  is  very  limited  due  to  frequency-dependent  attenua¬ 
tion;  second,  the  multi-path  delay  spread  is  excessively  long; 
third,  the  Doppler  spread  effect  is  significantly  large.  Due 
to  the  harsh  channel  condition,  demonstrated  rates  in  UWA 
communications  are  limited  to  several  kilo-bits  per  second 
within  40  km  range.  Multiple-input,  multiple-output  (MIMO) 
communications  which  improve  the  data  rate  significantly  by 
employing  parallel  independent  propagation  channels  existing 
in  the  ocean,  however,  are  much  harder  than  single-input, 
multiple-output  (SIMO)  communications  due  to  not  only  the 
inter-symbol  interference  (ISI)  but  also  the  multi-access  inter¬ 
ference  (MI)  from  signals  other  than  the  desired  signal. 

Many  approaches  have  been  proposed  for  UWA  re¬ 
ceiver  designs.  Earlier  UWA  systems  have  used  non-coherent 
frequency- shift  keying  (FSK).  Bandwidth  and  power  efficient 
coherent  modulation  scheme  has  been  first  developed  in  early 
1990’s  [1].  In  recent  years,  much  progress  has  been  made  [2]- 
[7]  in  coherent  UWA  communications.  The  frequency-domain 
(FD)  designs  using  single-carrier  frequency-domain  equaliza¬ 
tion  (SC-FDE)  and  orthogonal  frequency  division  multiplexing 
(OFDM),  are  found  in  [2]  and  [3],  [4],  respectively.  The 
time-domain  (TD)  receiver  designs  using  channel  equalization 
are  found  in  [5]  and  [6],  [7]  using  MIMO  technology.  The 
advantage  of  the  TD  designs  is  that  no  gaps  are  inserted  among 
transmission  blocks,  resulting  in  higher  data  efficiency  over 
common  FD  block  transmission  methods. 


The  advent  of  iterative  equalization  technology  [8]  has 
enabled  powerful  receiver  design.  The  iterative  equalization 
for  UWA  communications  has  been  first  proposed  in  [9]  with 
a  joint  maximum  a  posteriori  probability  (MAP)  iterative 
equalizer  and  turbo  decoder.  However,  the  complexity  of 
MAP  equalizer  is  prohibitively  high,  and  a  low-complexity 
iterative  decision-feedback  equalizer  (DFE)  has  been  proposed 
and  tested  by  a  multichannel  UWA  transmission  [10].  It¬ 
erative  DFE  for  MIMO  systems  has  been  proposed  in  [7], 
using  a  multi-band  transmission.  Iterative  equalization  using 
linear  equalizer  (LE)  [11]  has  also  been  applied  to  UWA 
communications  [12],  [13],  with  extensive  experimental  tests 
demonstrating  the  effectiveness  of  the  iterative  LE  scheme. 

Hard-decision  block  decision-feedback  equalizer  (BDFE), 
which  performs  better  in  fast  time-varying  multi-path  channel 
than  its  conventional  DFE  counterpart,  has  been  proposed 
in  [14],  and  soft-decision  iterative  BDFE  has  been  proposed 
in  [16],  both  for  single-input  single-output  (SISO)  systems. 
In  this  paper,  we  develop  an  iterative  BDFE  for  MIMO 
systems  and  apply  it  to  high-rate  UWA  communications. 
The  proposed  MIMO  BDFE  performs  successive  interference 
cancelation  (SIC)  of  both  ISI  in  the  time  domain  and  MI  in 
the  space  domain.  It  enables  a  near-optimum  sequence-based 
log-likelihood  ratio  (LLR)  evaluation,  which  is  superior  to  the 
conventional  symbol-based  LLR  calculations.  The  proposed 
receiver  design  is  tested  by  extensive  experimental  data  col¬ 
lected  at  the  Martha’s  Vineyard,  Edgartown,  MA,  in  October 
2008,  and  at  the  Gulf  of  Mexico,  in  July  2008. 

II.  System  Model 

Consider  a  MIMO  UWA  system  with  N  transducers  and  M 
hydrophones,  shown  in  Fig.  1,  where  an ^  denoting  the  k- th  bit 
of  the  n-th  stream,  is  encoded,  interleaved,  and  modulated,  into 
frn,/c,  cn?fc,  and  sn?fc.  The  modulated  symbols,  sn 5&,  are  format¬ 
ted  for  transmission.  For  the  system  employing  a  constellation 
of  size  Q,  the  group  of  coded  bits,  {cn,(k-i)\o&2Q+P}Xp=i  >  are 
mapped  to  the  modulation  symbol  sn^- 

The  baseband  signal  received  at  the  m-th  hydrophone  is 
expressed  by 

N  L  —  l 

Um^k  =  ^  ^  y  ^  ^)^n,k—l  H-  V m,k  (1) 

n=l 1=0 

where  hn,m(k,l)  is  the  l- th  fading  coefficient  between  the  n-th 
transducer  and  the  m-th  hydrophone  at  time  k. 
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Fig.  1.  MIMO  transmitter. 


III.  Iterative  MIMO  Receiver 


The  structure  of  the  proposed  iterative  MIMO  receiver  is 
depicted  in  Fig.  2,  where  the  soft-decision  MIMO  BDFE  iter¬ 
atively  exchanges  soft  decisions  with  MAP  decoders,  through 
the  interleavers  (n)  and  the  de-interleavers  (n-1)-  The  MIMO 
BDFE  is  performed  jointly  for  M  received  streams  to  yield 
N  soft-decision  streams  {A  (cn?/c)}^==1.  The  soft  decisions, 
with  the  a  priori  soft  decisions  being  subtracted  out,  are 
delivered  through  the  de-interleavers  for  MAP  decoding.  The 
MAP  decoders  generate  new  soft  decisions,  { AD  (6n,fe)}n_1, 
about  coded  bits.  The  new  soft  streams,  with  the  a  priori 
soft  decisions  being  subtracted  out,  are  fed  back  through 
the  interleavers  for  equalization.  Detection  performance  is 
gradually  improved  over  iterations  and  the  final  hard  decisions 
{a n,kjn=i  are  made  after  convergence. 


where  H  E  CMKxNK  is  a  block  matrix  with  the  (m,n)-th 
block  being  channel  matrix  Hm?n  defined  as 


'ftm,n(0)  0  0 


0 

0 


(3) 


0  ’  1)  ’ 


The  soft-decision  MIMO-BDFE  performs  detection  of  the 
k- th  symbol  from  the  n-th  data  stream,  sn^,  through  a 
feedforward  filter,  Wn?/c  E  CNKxMK ,  and  a  zero-diagonal 
upper  triangular  feedback  filter,  Bn?fc  E  CNKxNK,  as  shown 
in  Fig.  3.  The  equalized  symbol  vector  before  LLR  evaluation 


Fig.  3.  Soft-decision  MIMO-BDFE  (Wn5&  and  are  the  feedforward 

and  feedback  MIMO  equalizer  matrices  when  the  k- th  symbol  of  the  n-th 
stream  is  equalized.) 

can  be  written  as 


Sfi,k — ^^n,k  (y  Bn?/C  (s  Snj/C)  (4) 

where  the  estimated  channel  H  in  Fig.  3  is  replaced  by 
the  true  channel  H  for  the  convenience  of  denotation,  and 
sn  k  is  obtained  by  setting  the  mean  of  sn^  in  the  mean 
vector  s  =  E[s]  to  zero  while  leaving  all  the  other  ele¬ 
ments  unchanged.  The  operation  E[-]  denotes  mathematical 
expectation.  The  employment  of  sn^  instead  of  s  is  to  avoid 
the  instability  caused  by  positive  feedback.  The  vector  s  = 
[si,i,  -  *  -  ,SN,  l,-"  ,sN)KY  £  cNKxl  contains  the 

soft-decision  symbols. 

With  the  common  assumption  of  perfect  decision  feedback, 
i. e. ,  s  =  s,  the  error  vector  of  BDFE  can  be  written  as 


Fig.  2.  Iterative  receiver  using  soft-decision  MIMO  BDFE. 


A.  Soft-decision  MIMO-BDFE 

Consider  a  length-Ff  symbol  block,  sn  = 
Kir  -  ,«»,*]*  E  SKxl ,  where  (•)*  denotes  matrix 
transpose  and  S  =  {Xqi^Li  is  the  constellation  set  with 
cardinality  Q.  Within  each  block,  the  channel  is  treated 
as  time-invariant  [14].  Defining  received  signal  vector 
ym  =  [2/m,  1, 2/m, 2,  •  *  •  ,  Vm,K\,  additive  noise  vector 

for  the  ra-th  hydrophone  and 
stacking  up  {ym}^=1,  {vm}^=1  and  {sn};)=1  into  column 
vectors  leads  to 


en,/c=Wn?/c  (y-Hsn?/c)-(Bn?/c+lATK)  (s-sn^)  (5) 

where  Ink  is  an  identity  matrix  of  size  NK.  Minimizing  the 


mean  square  error,  E 


feedforward  and  feedback  matrices 


en,ken,k 


,  leads  to  the  solutions  for  the 


B n,k  —  Gn,k  ~  W,  (6a) 

W„>fc  =  Gn,k$n,kHh  [H$n>fcH'1  +  <j2vInk]  _1  (6b) 

where  (-)h  denotes  matrix  conjugate  transpose,  <&n,k  is  ob¬ 
tained  by  setting  the  variance  of  sn^  in  the  covariance  matrix, 
=  E  [(s  —  s)(s  —  s)h]  as  one,  and  G n,k  C  CNKxNK  is 
an  unit-diagonal  upper  triangular  matrix  obtained  from  the 
Cholesky  decomposition  of  +  ^2  H^H  as 

+  4iFH  =  G^Dn,kG„,t 

Oy 


y  =  Hs  +  v 


(2) 


(7) 


with  D Ujk  G  CNKxNK  being  a  diagonal  matrix.  With  the 
filters  in  (6),  the  original  system  is  converted  to  an  equivalent 
system  as  (c.f.  (5)) 

rn,fc  -  Wn?,  (y  -  Hsnjfc)  =  Gn?/c  (s  -  sn>fe)  +  enjfc.  (8) 


B.  Sequence-based  LLR  Calculation 

In  the  i- th  iteration,  the  sequence-based  a  posteriori  prob- 

(i) 

ability  (APP)  of  snJk  conditioned  on  rn ^  is 


(Tn,k\8Sk)  P  (S«,fc) 

p(rn,k) 


(9) 


Based  on  the  assumption  that  en ^  in  (8)  is  a  random 
vector  with  independent  Gaussian  elements,  we  can  simplify 
the  probability  calculation  as 


P 


NK 

)  ~  TT  __Pxn  < 

Pn,k 

2 

j  ~  1  1  exP  S 

/  ±±  TTCFj 

3  =  1 

f 

-I 

_ 

(10) 


where  cr2  =  djj  with  djj  being  the  j -th  diagonal  element 
of  D n,k>  and  PJn  k  defined  in  (11)  at  the  top  of  the  next 
page.  In  (11),  gjj  is  the  (j,  Z)-th  element  of  Gn^,  rj  is  the 
j-th  element  of  rn^,  J  =  (fc  — 1)7V+  n ,  and  the  operations 
p(l)  =  \l/K ]  and  q(l)  =  l  —  (\l/K]  —  l)K  represent  the  mapping 
from  the  linear  index  set,  {l}iK,  to  the  double  index  set, 
The  tentative  soft  decisions  sp^  and  s$q  are 
taken  from  the  (i  —  l)-th  iteration  and  the  i-th  iteration, 
respectively.  The  upper  triangular  structure  of  Gn^  determines 
that  the  symbols,  sp,q,  with  index  satisfying  (p  —  1  )K  +  q  > 
(n  —  1)K  +  k,  are  detected  before  sn^,  thus  the  tentative 
detections  of  those  symbols  from  the  current  iteration  can  be 
used.  The  probabilities  P  (s^k)  and  P(rn,k)  are  determined 
by  the  a  priori  probabilities  or  mapping  bits,  and  the  nor¬ 
malization  Y^=ip{^n,k  =  X<zlrn,/c)  =  1,  respectively.  With 

P  the  soft-decision  symbols  can  be  calculated, 

and  the  APP  and  LLR  of  the  mapping  bits,  are  evaluated  as 
follows 


P 


'  W 

Cn,(/c-l)log2<2+p 


b\rntk 


E  P{SnPn,k)  (12) 

S^kes{pb) 


where  Sp ^  =  {Xq\Xq  ^  <5  :  cq#  =  b}  with  b  G  {0, 1},  and 
p  =  1,  •  •  •  ,  log2  Q.  The  a  priori  LLR  is  subtracted  from  (13) 
to  produce  extrinsic  LLR  for  MAP  decoders. 


C.  Low -complexity  Approximate  Soft-decision  MIMO-BDFE 

The  major  computational  burden  of  the  proposed  iterative 
equalization  comes  from  the  requirement  of  updating  the  filter 
matrices  W n,k  and  &n,k  for  each  symbol.  In  each  updating, 
two  matrix  inversions  and  one  Cholesky  decomposition  are 
involved.  A  closer  observation  reveals  that  the  symbol-wise 
filter  updating  is  solely  due  to  the  dependence  on  the  second- 
order  a  priori  information  3>n?/c.  Therefore,  the  computational 
complexity  can  be  reduced  significantly  by  replacing  &n,k 
with  <i>,  which  leads  to  constant  filter  matrices,  W  and  B, 
for  all  symbols  in  the  block.  It  is  verified  by  both  simulation 
and  data  processing  that  employing  constant  filter  matrices 
for  all  symbols  doesn’t  apparently  degrade  the  equalization 
performance. 

The  proposed  sequence-based  equalization  gleans  bit  relia¬ 
bility  information  from  the  whole  sequence  rn^  (or  equiv¬ 
alently,  s n,k),  rather  than  the  single  symbol  sn^k  used  by 
the  symbol-based  iterative  equalization,  it  thus  achieves  better 
performance.  The  sequence-based  LLR  calculation  is  more 
computationally  intensive.  However,  the  evaluation  can  be 
simplified  by  utilizing  the  sparse  structure  of  the  equivalent 
channel  matrix  Gn?fc.  Details  are  referred  to  [16]  for  SISO  sys¬ 
tems,  and  similar  results  apply  to  MIMO  systems.  Generally, 
the  proposed  equalization  has  the  same  order  of  computational 
complexity  as  the  symbol-based  iterative  equalizations  in  each 
iteration,  while  gains  more  overall  computational  savings  due 
to  its  faster  converges. 

IV.  Experimental  Results 

The  proposed  iterative  receiver  using  MIMO  BDFE  was 
tested  by  two  undersea  trials  of  UWA  communications.  The 
first  trial  named  SPACE08,  was  conducted  at  the  coast  of 
Martha’s  Vineyard,  Edgartown,  MA,  in  October  2008.  The 
second  trial  named  GOMEX08,  was  conducted  at  the  Gulf  of 
Mexico,  in  July  2008.  Details  on  the  detections  of  real  data 
are  presented  in  the  following  two  subsections. 

A.  Results  of  SPACE08  experiment 

In  this  experiment,  the  symbol  rate  was  9.765625  kilo- 
symbols  per  second.  The  carrier  frequency  was  13  kHz. 
Modulation  schemes  included  QPSK,  8PSK  and  16QAM.  The 
water  depth  was  about  15  meters.  The  transmitting  equipment 
consisting  of  four  transducers  located  about  4  meters  above 
the  sea  bottom.  The  receiving  equipment  consisted  of  six 
hydrophone  arrays,  denoted  as  SI  (60  m,  Southeast),  S2  (60 
m,  Southwest),  S3  (200  m,  Southeast),  S4  (200  m,  Southwest), 
S5  (1  km,  Southeast)  and  S6  (1  km,  Southwest),  respectively. 
All  received  files  contain  12  channels  of  data  named  after  their 
corresponding  receiving  arrays. 

Fig.  4  depicts  one  transmission  frame,  which  contains  the 
head  linear  frequency  modulation  (LFM)  signal  named  LFMB, 
three  packets  with  QPSK,  8PSK  and  16QAM  modulations,  and 
the  tail  LFM  signal  named  LFME.  The  LFMB  and  LFME 
signals  aim  for  the  frame  (and  packet)  coarse  synchronization 
and  the  channel  length  measurement.  Each  packet  is  made  up 
of  a  m-sequence  of  length  511,  and  the  data  payload  consists 
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Fig.  4.  Frame  structure. 


of  30,000  modulation  symbols.  The  ra-sequence  plays  similar 
role  as  LFM  signal.  The  data  payload  consists  of  both  training 
blocks  and  information  blocks,  with  no  gaps  among  blocks. 
The  block  size  is  designated  during  detection,  depending  on 
the  channel  condition.  The  placement  of  training  blocks  is  also 
flexible,  depending  on  the  tradeoff  during  detection. 

The  iterative  receiver  applied  to  experimental  packet  de¬ 
tection  is  described  as  follows.  As  mentioned  above,  the 
received  signal  was  artificially  partitioned  into  blocks  with 
small  time  durations.  We  chose  the  block  size  as  =  200. 
The  MIMO  channel  with  measured  channel  length  of  L  =  100, 
was  estimated  either  in  training  mode  or  in  decision-directed 
(DD)  mode.  In  training  mode,  the  channel  was  estimated 
with  three  training  blocks  (or  Np  =  600  training  symbols), 
inserted  in  the  data  payload  for  every  T  blocks,  where  T 
was  properly  selected  so  that  tradeoff  between  the  detection 
performance  and  the  training  overhead  was  achieved.  In  DD 
mode,  Np  =  600  previously  detected  symbols  are  used  to  re- 
estimate  the  channel,  for  equalizing  the  current  block.  To  apply 
BDFE  presented  in  Section  III-A,  the  post-cursor  effect  from 
previous  block  was  removed  from  the  current  block,  using  the 
equalized  soft-decision  symbols  and  the  estimated  channel.  To 
account  for  the  signal-to-noise  ratio  (SNR)  degradation  in  the 
equalized  symbols  at  the  tail  of  current  block  [14],  overlapped 
block  detection  was  used,  for  which  some  equalized  symbols  at 
the  tail  of  current  block  was  detected  again  in  the  next  block. 
During  the  detection,  we  selected  an  overlapping  length  of 
TVovip  =  50  which  was  conservative  since  shorter  overlapping 
length  still  worked  well.  The  percentage  of  the  training  over¬ 
head  could  be  roughly  estimated  as  r\  —  T(N^-n0Viv)  x  100%. 

Fig.  5  shows  an  example  of  the  estimated  channel  impulse 
responses  (CIRs)  for  the  200  m  two-transducer  QPSK  trans¬ 
missions.  It  is  obvious  from  the  figure  that  the  subchannels 
are  sparse  and  some  are  non-minimum  phase,  which  imposes 
difficulty  on  robust  UWA  communication. 

The  scatter  plots  in  Fig.  6  and  Fig.  7  show  the  soft-decision 
symbols  at  the  outputs  of  the  equalizer  and  the  MAP  decoder 
for  QPSK  and  8PSK,  respectively.  The  improvement  over 
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Fig.  5.  Estimated  channel  impulse  responses  for  QPSK  transmission. 


iterations  is  clearly  observed  by  comparing  subfigures  in  each 
row  of  both  figures.  The  performance  improvement  after  MAP 
decoder  is  also  observed  by  comparing  subfigues  in  each 
column  of  both  figures. 


Fig.  6.  Soft-decision  symbols  at  the  outputs  of  the  soft-decision  equalizer 
and  the  MAP  decoder  for  QPSK  modulation. 


For  200  m  transmission,  45  packets  (30  S3  packets  and  15 
S4  packets)  were  processed.  The  results  for  the  two-transducer 
QPSK  transmission  are  listed  in  Table  I,  where  T  =  30 
is  selected,  incurring  the  training  overhead  of  14%.  All  45 
packets  achieve  zero  bit  error  rates  (BERs),  for  which  most 


TABLE  III 


BER  OF  3  X  12  MIMO  1000  M  QPSK  TRANSMISSION 


Transducer 

Iter.  1 

Iter.  2 

Iter.  3 

Iter.  4 

1 

9.333e-4 

0 

0 

0 

2 

3.053e-2 

2.333e-4 

0 

0 

3 

8.927e-2 

1.800e-3 

0 

0 

TABLE  IV 


BER  OF  4  X  12  MIMO  1000  M  QPSK  TRANSMISSION 


Transducer 

Iter.  1 

Iter.  2 

Iter.  3 

Iter.  4 

1 

1.318e-3 

0 

0 

0 

2 

1.600e-2 

1.364e-4 

0 

0 

3 

6.577e-2 

1.455e-3 

0 

0 

4 

1.100e-2 

0 

0 

0 

Fig.  7.  Soft-decision  symbols  at  the  outputs  of  the  soft-decision  equalizer 
and  the  MAP  decoder  for  8PSK  modulation. 


TABLE  I 

Results  for  2  x  12  MIMO  200  m  QPSK  transmission 


Number  of  iterations 
to  achieve  zero  BER 

Number  of 
packets 

1 

20 

2 

15 

3 

6 

4 

4 

packets  only  require  two  iterations.  For  the  two-transducer 
8PSK  transmission,  25  out  of  45  packets  achieve  zero  BERs 
within  five  iterations.  The  average  BERs  for  remaining  20 
packets  are  5.657  x  10-4  and  6.018  x  10-4,  respectively,  for 
the  first  transducer  and  the  second  transducer.  The  results  for 

the  1000  m,  two-transducer,  QPSK  transmission  are  listed  in 
Table  II,  where  T  =  20  incurring  20%  training  overhead,  is 
selected.  All  15  processed  packets  (5  S5  packets  and  10  S6 
packets),  achieve  zero  BER  after  three  iterations.  For  8PSK 
case,  we  get  5  packets  achieving  zero  BERs  within  three 
iterations,  and  the  average  BER  for  the  remaining  10  packets 
are  8.024  x  10-5  and  1.191  x  10-3  for  the  first  transducer  and 
the  second  transducer,  respectively. 

We  next  demonstrate  the  results  for  the  1000  m  three- 
transducer  and  four-transducer  QPSK  transmissions.  The  BER 
results  are  listed  in  Table  III  and  Table  IV,  with  the  training 
overhead  being  24%  and  34.7%,  respectively.  For  both  cases, 
zero  BERs  are  achieved  after  three  iterations. 


B.  Results  of  GOMEX08  experiment 

In  this  experiment,  a  moving  transmitting  equipment  and 
a  fixed  receiving  equipment  was  deployed.  The  transmitting 
equipment  consisted  of  four  transducers  located  about  50  m 
deep  in  the  sea  at  zero  speed,  and  rose  to  30  m  at  high 
speed  of  2.7  knots.  The  receiving  equipment  consisted  of  eight 
hydrophones  was  fixed  50  m  deep  in  the  sea.  The  transmission 
distance  was  1.7  km  to  2  km  and  the  symbol  rate  was  4  kilo- 
symbols  per  second.  The  carrier  frequency  was  17  kHz.  Each 
transmission  frame  consisted  of  multiple  packets.  The  number 
of  packets  within  one  frame  was  48,  24  and  12  corresponding 
to  the  different  packet  sizes  of  1024,  2048  and  4096.  The 
channel  length  was  measured  as  L  =  80  in  this  trial,  and  a 
block  size  of  N \  =  150  was  adopted  during  detection.  The 
overlapping  length  was  selected  as  7Vovip  =  40.  An  example 
of  the  estimated  channel  is  demonstrated  in  Fig.  8,  where  the 
sparseness  and  non-minimum  phase  of  the  sub-channels  are 
observed,  similar  to  Fig.  5. 


CIR:  Tx1  to  Rx1  CIR:  Tx1  to  Rxg 
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TABLE  II 

Results  for  2  x  12  MIMO  1000  m  QPSK  transmission 


Number  of  iterations 
to  achieve  zero  BER 

Number  of 
packets 

1 

9 

2 

5 

3 

1 

Fig.  8.  Estimated  channel  impulse  responses  for  GOMEX08  experiment. 

We  present  the  detection  results  for  three  frames  with  packet 
size  2048  and  three  frames  with  packet  size  4096,  both  with 
QPSK  modulations.  The  number  of  iterations  required  for 
achieving  zero-BER  for  each  packet,  are  listed  in  Table  V  and 


TABLE  V 

Number  of  iterations  for  achieving  zero  BER:  2x8  MIMO 
QPSK  TRANSMISSION  WITH  PACKET  SIZE  2048 


Index  of  packet 

Frame  1 

Frame  2 

Frame  3 

1 

1 

1 

1 

2 

2 

1 

2 

3 

1 

1 

1 

4 

1 

1 

1 

5 

1 

1 

1 

6 

1 

1 

1 

7 

2 

1 

1 

8 

1 

1 

1 

9 

1 

1 

2 

10 

1 

1 

3 

11 

1 

3 

1 

12 

1 

2 

1 

13 

1 

1 

1 

14 

1 

1 

1 

15 

1 

1 

3 

16 

1 

1 

3 

17 

1 

1 

1 

18 

1 

4 

1 

19 

1 

1 

1 

20 

1 

1 

1 

21 

2 

1 

1 

22 

1 

1 

1 

23 

1 

1 

1 

24 

1 

2 

1 

TABLE  VI 

Number  of  iterations  for  achieving  zero  BER:  2x8  MIMO 
QPSK  TRANSMISSION  WITH  PACKET  SIZE  4096 


Index  of  packet 

Frame  1 

Frame  2 

Frame  3 

1 

4 

1 

1 

2 

2 

2 

1 

3 

1 

3 

1 

4 

1 

2 

1 

5 

4 

2 

2 

6 

3 

1 

5 

7 

1 

1 

3 

8 

4 

1 

2 

9 

2 

3 

2 

10 

1 

3 

1 

11 

1 

1 

2 

12 

3 

2 

1 

Table  VI.  From  both  tables,  most  of  the  packets  require  no 
more  than  three  iterations  to  achieve  zero  BER,  which  again 
verify  the  robustness  and  fast  convergence  of  the  proposed 
turbo  receiver. 


V.  Conclusion 

We  have  demonstrated  a  new  receiver  structure  for  high- 
rate  UWA  communications,  using  iterative  MIMO  BDFE.  The 
proposed  soft-decision  MIMO  BDFE  has  performed  succes¬ 
sive  interference  cancelation  of  both  ISI  in  the  time  domain 
and  MI  in  the  space  domain,  and  it  has  enabled  a  sequence- 
based  LLR  calculation  leading  to  excellent  detection  perfor¬ 
mance  with  fast  convergence.  The  MIMO  channel  has  been 
estimated  and  tracked  relying  on  either  training  symbols  or 


previously  detected  symbols.  Rigorous  testing  of  the  proposed 
receiver  scheme  by  extensive  experimental  data  collected  in 
the  SPACE08  and  GOMEX08  undersea  trials,  has  shown 
that  consistent  reliable  detection  was  achieved  with  moderate 
training  overhead. 
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